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ABSTRACT 

We report the detection of 158 /xm [CII] fine-structure line emission from MIPS J142824. 0-1-352619, 
a hyperluminous (Ljr ~ 10^'^ Lq) starburst galaxy at ^ = 1.3. The line is bright, and corresponds to 
a fraction L[cii]/iFiR ~ 2 x 10"'^ of the far-IR (FIR) continuum. The [CII], CO, and FIR continuum 
emission may be modeled as arising from photodissociation regions (PDRs) that have a characteristic 
gas density of n ~ lO'*'^ cm"'^, and that are illuminated by a far-UV radiation field ~10'^'^ times more 
intense than the local interstellar radiation field. The mass in these PDRs accounts for approximately 
half of the molecular gas mass in this galaxy. The L[cjj]/Lfir ratio is higher than observed in local 
ULIRGs or in the few high-redshift QSOs detected in [CII], but the L[cii]/^fir and ico/^FiR ratios 
are similar to the values seen in nearby starburst galaxies. This suggests that MIPS J142824. 0-1-352619 
is a scaled-up version of a starburst nucleus, with the burst extended over several kiloparsecs. 
Subject headings: galaxies: individual(MIPS J142824. 0+352619) — galaxies: high-redshift — galaxies: 
ISM — galaxies: starburst — infrared: galaxies — submillimeter 



1. INTRODUCTION 

The 2P3/2 ^Pi/2 transition of C+ (A = 157.74 /im) is 
one of the brightest emission lines in star-forming galax- 
ies, typically accounting for 0.1 — 1% of the far-IR (FIR) 
continuum (Stacey et al. 1991; Malhotra et al. 2001). 
Much of this emission arises from the warm and dense 
photodissociation regions (PDRs) that form on the UV- 
illuminated surfaces of molecular clouds. The [CII] tran- 
sition is a primary PDR coolant, and is a sensitive probe 
of both the physical conditions of the photodissociated 
gas, and the intensity of the ambient stellar radiation 
field (Hollenbach & Tielens 1999). In ultraluminous in- 
frared galaxies (ULIRGs), the L[cii]/Lfir, ratio is a fac- 
tor of ^7 times lower than in less luminous systems (Luh- 
man et al. 2003), for reasons that are not well under- 
stood. 

Large aperture (sub)millimeter telescopes have been 
used to search for [CIl] emission from sources at z > 3, 
and 3 FIR-luminous QSOs at z = 4.4 — 6.4 have been 
detected thus far. SDSS J1148 (Maiohno et al. 2005) and 
BR 1202N (lono et al. 2006c) have low L[cii]/Lfir ratios 
similar to or smaller than the mean local ULIRG value, 
while BRI 0952 (Maiolino et al. 2009) has a somewhat 
larger ratio falling at the lower end of the range seen 
in normal galaxies. We have initiated a search for [CII] 
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emission from galaxies at z — 1—2, concentrating on star- 
formation-dominated systems selected by their FIR or 
submillimeter continuum brightness. Here we report our 
first detection, from the z = 1.3 hyperluminous starburst 
galaxy MIPS J142824.0+352619 (hereafter MIPS J1428). 

MIPS J1428 was identified in a Spitzer/MIPS blank 
field survey as a bright 160 /zm source with red 
optical/near-IR colors, and subsequent spectroscopic and 
photometric measurements established it as a hyperlu- 
minous (Lir[8 - 1000 Mm] = 3.2 x lO^^ ^^-^ galaxy at 
z — 1.325 (Borys et al. 2006). There are several indica- 
tions that this IR emission is powered by star formation, 
including the brightness of the PAH features and the non- 
detection in hard X-ray emission, and with a starburst 
origin the estimated Lir, corresponds to a star forma- 
tion rate of -5500 Mq yr'^ (Borys et al. 2006; Desai 
et al. 2006). MIPS J1428 was detected in C0(2^1) and 
C0(3— >2), and the line luminosities correspond to a large 
gas mass of Mh^ 10" Mq (lono et al. 2006b). How- 
ever, a gravitational lens may amplify the flux by as much 
as a factor of 10 (Borys et al. 2006; lono et al. 2006a, b), 
making the actual star formation rate and gas mass cor- 
respondingly lower. Independent of the unknown lensing 
amplification, the large value of im/i'co ~ 320 Lq (K 
km s~i)"i (lono et al. 2006b) indicates that MIPS J1428 
is forming stars with high efficiency, similar to that seen 
in local ULIRGs and high-redshift submillimeter galaxies 
(SMGs) (Solomon & Vanden Bout 2005). 

This is the first detection of [CII] in the z = 1 — 3 epoch 
of peak star formation activity in the Universe, and the 
first detection from a high-redshift galaxy not associated 
with a QSO. 



2. OBSERVATIONS 

We observed the [CII] line toward MIPS J1428 in 
March 2008 with ZEUS (Stacey et al. 2007; Hailey- 
Dunsheath 2009) at the Caltech Submillimeter Obser- 
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Spectral Line and Continuum Observations of MIPS J1428 
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Fig. 1. — ZEUS/CSO spectrum of the 158 ^m [CII] transition in 
MIPS J1428. The velocity scale is referenced to z = 1.325, and the 
pixel centered aXv = 1250 km s^^ is removed due to the presence of 
a telluric absorption feature. The integrated intensity is obtained 
by summing the 2 central spectral bins, which cover Av = 416 km 
. The line is detected at &a significance. 

vatory (CSO)** on Mauna Kca, Hawaii. ZEUS is a dual- 
band (350 and 450 fim) grating spectrometer equipped 
with a 1 X 32 semiconductor bolometer array. For this 
observing run the instrument was configured as described 
in Hailey-Dunsheath et al. (2008). At the redshifted 
wavelength (A ~ 366.75 /im) of the [CH] transition the 
spectrometer provides a slit-limited resolving power of 
A/ A A ~ 900, and the 14 detectors operating in the 350 
fjbva. window combine to provide an instantaneous band- 
width corresponding to 5v = 2825 km s~^. 

Orion (BN-KL) was used as spectral calibrator, and 
Mars (Tb = 213 K; Hildebrand et al. 1985) as flux 
calibrator. We integrated on MIPS J1428 for a total 
of 135 minutes (including chopping) in good submil- 
limeter weather, with Tsso^jm = 1-0 — 1.5. The 13" 
ZEUS/CSO beam was centered at R.A. = 14'^28'"24!1, 
decl. = +35°26'19" (J2000.0), and the data were ob- 
tained by chopping and nodding the telescope with a 
30" throw. A linear baseline is removed from the final 
spectrum, shown in Figure 1. The integrated line flux is 
F^Av = 713 Jy km s"'^, and we compare this with other 
observations in Table 1. With Lfir = 2.6 x 10^^ Lq (Bo- 
rys et al. 2006), we estimate Ljchj/Lfir ~ 2.1 x 10~^. 

3. RESULTS 

3.1. Origins of the [CII] Line 

Carbon has a lower ionization potential (11.26 eV) 
than hydrogen, and consequently the [CII] transition is 
an important coolant of both ionized gas and difi^use 
atomic gas. However, most of the [CII] emission from 
nearby IR-bright galaxies arises from dense PDRs il- 
luminated by the far-UV (FUV) radiation from young 
stars (Crawford ct al. 1985; Stacey et al. 1991). Com- 
bined HII region and FDR modeling of the starburst tem- 
plates NGC 253 (Carral et al. 1994) and M82 (Lord et al. 
1996; Colbert et al. 1999) has demonstrated that PDRs 
account for at least 70% of the [CII] emission in these 
sources. The available data suggests that such starburst 
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following Ffir = 1.26 x 10"" (2.58/60 -t- /lOo) W m'^ (Sanders 
& Mirabel 1996). 

galaxies arc the best local analogs of MIPS J1428 (see 
section 4) , so here we assume that 70% of the [CII] emis- 
sion from MIPS .11428 arises in PDRs. We fiirther as- 
sume that PDRs produce afl of the C0(2^1), CO(3^2), 
and FIR continuum emission in this source, and we model 
these tracers as arising from a single representative PDR. 

3.2. PDR Analysis 

We use the PDR models of Kaufman et al. (1999), 
in which the emission is determined by the gas density, 
n, and the incident FUV (6 cY < hu < 13.6 eV) flux, 
Go (expressed in units of the Habing Field, 1.6 x 10~^ 
ergs cm~^ s~^). These models calculate the line emission 
generated by a cloud illuminated on only one side, and 
must be modified to model clouds in the active regions 
of galaxies that are illuminated on all sides. For this ge- 
ometry, Kaufman ct al. (1999) note that an observer will 
detect optically thin radiation emitted by both the near 
and far sides of the cloud, while optically thick emission 
will only be visible from the near side. The [CII] line is 
generally moderately optically thin in these calculations, 
and the FIR continuum is also assumed to be optically 
thin, but the low- J CO transitions are optically thick for 
the nominal Ay = 10 cloud depth. We therefore increase 
the measured CO fluxes by a factor of 2 with respect to 
these other tracers before comparing with the models. 

In Figure 2 we show the values of n and Go allowed 
by the observations listed in Table 1, after corrcicting the 
[CII] and CO measurements as described above. In the 
n « 10^ — 10^ cm~^ range the L^qh]/ L-piji ratio provides 
a strong constraint on Gq, while the i[cii]/-^co(2^-i) 
ratio is primarily a function of density. The best fit 
solution obtained from these two ratios is consistent 
with the n > lO"^'^ cm~^ lower limit imposed by the 
^co(3->2)/-^co(2^i) ratio, and we conclude that the rel- 
ative strengths of the [CII], C0(2^1), CO(3-S>2), and 
FIR continuum emission may be reproduced by a PDR 
with n ~ 10^-2 cm-3 and Go ~ 10^-^. Reducing the PDR 
contribution to the total [CII] emission by a factor of 2 
(3) would increase n by a factor of 7.2 (22), but Go by 
only a factor of 1.6 (1.8). We also note that without the 
factor of 2 correction to the observed CO fluxes, n would 
be 4.2 times lower, but Go would remain unchanged. 

An estimate of the atomic gas mass associated with 
the PDRs in MIPS J1428 is obtained from the inferred 
[CII] luminosity and this PDR solution. Assuming the 
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Fig. 2. — i[cii]/-f'FiR, ■f'[Cii]/-^^CO(2^i)> and 

^CO{3->2)/icO(2^i) as a function of density (n) and FUV 
flux (Go), from the PDR models of Kaufman et al. (1999). 
Contours are drawn for the measured values {solid) and Ha 
uncertainties {dot-dashed) of MIPS J1428, after scaling the [CII] 
flux by 70% to account for a non-PDR contribution, and increasing 
the CO fluxes by a factor of 2 to account for geometrical effects 
(see section 3). 

[CII] emission is optically thin, this mass is given as 

/ 1 + 2 exp(-91 K/T) + n,,it/n \ 

2exp(-91K/T) j' ^' 

where we use a C"*" abundance per hydrogen atom of 
Xc+ = 1.4 X 10-4 (Savage & Sembach 1996). The sur- 
face temperature of a PDR with n — lO'*'^ cm"'^ and 
Go = 103-2 is Ts « 230 K (Kaufman et al. 1999). As- 
suming this temperature is representative of the full C+ 
region, and adopting a critical density of ricrit = 2.7 x 10"^ 
cm~^ (Launay & RouefF 1977), this expression yields 
Ma = 5.5 X 10^"/^-^ Mq (where fj, is the lensing am- 
plification). This is 55% of the molecular gas mass esti- 
mated from the CO luminosity (lono et al. 2006b) , which 
is similar to the PDR mass fractions in starburst galax- 
ies (Stacey et al. 1991, corrected for the smaller Xq+ 
used here). Much like in nearby starburst nuclei, the 
PDR mass in MIPS J1428 is a significant fraction of the 
total molecular gas mass. 

4. DISCUSSION 

4.1. [CII], CO, and the FIR Continuum: A 
Comparison to Other Calaxies 

To place our [CII] measurement of MIPS J 1428 in con- 
text, we compare the [CII], CO, and FIR continuum 
emission from this source with that from other galaxies 
and Galactic sources. In Figure 3 we plot the L[cii]/iFiR 
and icO(i->o)/^FiR ratios for samples of Galactic star- 
forming regions and nearby gas-rich galaxies (Stacey 
et al. 1991), normal galaxies (Malhotra et al. 2001), local 
ULIRGs (Luhman et al. 2003), and MIPS J1428 and the 
3 other high-rcdshift sources detected in [CII] (Maiolino 
et al. 2005; lono et al. 2006c; Maiohno et al. 2009). We 
estimate the CO(1->0) luminosity of MIPS J1428 assum- 
ing an i'co(2-!-i)/-^'co(i-!-o) = 0.9 ratio intermediate be- 
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tween the ratios observed in local ULIRGs (~0.7; Downes 
& Solomon 1998) and in starburst nuclei («1.1; Harrison 
et al. 1999; Weifi et al. 2005). For reference, we overplot 
the PDR model curves from Kaufman et al. (1999). A 
vector attached to the MIPS J1428 data point indicates 
how the measured ratios were adjusted for the PDR anal- 
ysis in section 3.2 (to account for non-PDR contributions 
to the [CII] emission, and optical depth effects in the CO 
transititions), and we further note that MIPS J1428 is 
placed in Figure 3 assuming an i'co(2->i)/^'co(i->o) i'^- 
tio ~20% lower than our best fit PDR solution. We com- 
pare MIPS J1428 with each of the populations shown in 
Figure 3 in turn. 

Stacey et al. (1991) study the [CII] emission from 
the central regions of a sample of nearby galaxies, and 
find that the i[cn]/-^^co(i->o) and i[cii]/iFiR ratios are 
sensitive tracers of the star formation activity in these 
sources. Starburst galaxies are well characterized by the 
same i[cii]/^co(i->o) = 4100 ratio seen in Galactic star- 
forming regions, which is a factor of ^3 times larger than 
in non-starburst galaxies and Galactic molecular clouds 
(Fig. 3). However, the i[cn]/iFiR ratios in starburst 
galaxies are higher than those observed in the most in- 
tense Galactic OB star formation regions. These results 
imply that while much of the molecular gas in a starburst 
nucleus is photodissociated by the FUV radiation from 
young stars, the very intense fields found in the most 
extreme Galactic star-forming regions are not providing 
the bulk of the luminosity. This picture is supported 
by PDR modeling, which indicates that the FUV fields 
producing the emission in the centers of starburst galax- 
ies are intermediate in intensity between those producing 
the emission in non-starburst galaxies, and those found 
in Galactic OB star formation regions (Fig. 3). As the 
data point for MIPS J1428 faUs near the high-Go end of 
the distribution of starburst points in Figure 3, we con- 
clude that in MIPS J1428, like in nearby starburst nuclei, 
the bulk of the molecular gas is exposed to moderately- 
intense FUV radiation produced by young stars. 

Malhotra et al. (2001) observe [CII] and other FIR fine- 
structure line emission from a sample of 60 normal galax- 
ies with varying levels of star formation activity, which 
are all sufficiently distant that the FIR emission is con- 
tained within the ?»70" ZSO-LWS beam. These authors 
find that the FIR continuum and much of the [CII] emis- 
sion arises from PDRs, which in many of their galaxies 
are characterized by similarly elevated values of Go as de- 
duced here for MIPS J1428. However, the Malhotra et al. 
(2001) sources have lower i[cii]/ico ratios than MIPS 
J1428 or its starburst analogs, indicating that most of the 
CO emission is produced by molecular gas residing in less 
active star-forming regions than those that produce the 
fine-structure line emission. This difference helps rein- 
force our conclusion that in MIPS J1428 it is the entire 
galaxy, rather than just an active central region, that is 
host to vigorous star formation. 

The median Ljcuj/LpiR ratio in local ULIRGs is a fac- 
tor of ^7 times lower than in normal star-forming galax- 
ies (Luhman et al. 2003), and a factor of '^4 times lower 
than in MIPS J1428. One possible explanation for this 
low ratio is that the [CII] emission is produced in dense 
PDRs illuminated by intense FUV radiation (e.g., Pa- 
padopoulos et al. 2007), as indicated by the PDR model 
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Fig. 3. — L[cjj]/LpiR vs. icO(i— >0)/-^FlR Galactic star-forming regions (crosses), starburst nuclei (filled squares), non-starburst 
nuclei (open squares), normal galaxies (triangles) , local ULIRGs (circles), and high-redshift sources (asterisks) . CO(l— >-0) luminosities of 
the normal galaxies and local ULIRGs are taken from the literature and from J. Gracia-Carpio et al. (2010, in preparation). CO(l— >0) 
luminosities for MIPS J1428 and 2 of the other high-redshift sources are estimated from the measurements of higher- J lines as described 
in the text. Starburst galaxies and Galactic star-forming regions are characterized by i[cil]/^CO(l^0) ~ 4100 (dashed). Overplotted are 
the PDR model calculations for gas density (n) and FUV field strength (Go) from Kaufman et al. (1999), and a vector indicates how the 
MIPS J1428 data point was shifted for the PDR analysis in section 3.2. 



overlays in Figure 3. Alternatively, the [CII] emission 
may originate in PDRs with more modest values of n and 
Go J but an additional source of FIR continuum emission 
may lower the global L[cii]/iFiR ratio (Luhman et al. 
2003; Abel et al. 2009). Regardless of the origins of the 
weak [CII] emission in local ULIRGs, it is clear that these 
galaxies provide poor templates for MIPS J1428. 

We also include the 3 high-redshift QSOs detected 
in [CII] in Figure 3. SDSS J1148 has been detected 
in CO (3-^2) and BRI 0952 in CO (5^4), and we use 



the measured i'co(3-f2)/^' 
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CO(l-fO) 



= 0.72 and modeled 



= 0.58 ratios for Mrk 231 (Pa- 



padopoulos et al. 2007) to estimate the C0(1— >0) lumi- 
nosities for these sources. To account for the uncertain- 
ties in the CO excitation of high-redshift systems (e.g., 
Hainline et al. 2006), we show these data points with 
errors bars corresponding to 0.2 dex and 0.3 dex, re- 
spectively. The low i[cii]/iFiR and ico/^FiR ratios in 
SDSS J1148 and BR 1202N indicate similar physical con- 
ditions as in local ULIRGs. Indeed, Maiolino et al. (2005) 
found that the [CII], CO, and FIR continuum emission 
of SDSS J1148 is consistent with a high-density, high- 
Go PDR, and Walter et al. (2009) showed that the star 
formation in this system is confined to a compact re- 
gion with a high surface brightness comparable to that 



at the center of Arp 220. However, BRI 0952 has a larger 
Lpiij/ipiR ratio and a i[cii]/-^^CO ratio similar to that 
in MIPS J1428. We suggest that while SDSS J1148 and 
BR 1202N are hkely dissimilar to MIPS J1428, the con- 
nection drawn here between MIPS J1428 and local star- 
burst nuclei may also provide insight into the nature of 
the host galaxy of BRI 0952. 

Finally, it is important to note that while we assume 
that the L^hj/Lfir and Lco/^fir ratios in MIPS J1428 
are unaffected by the potential gravitational lensing, we 
cannot rule out the possibility that these ratios are al- 
tered by differential magnification. For example, one 
could construct a model in which the intrinsic emis- 
sion ratios of MIPS J1428 are similar to those in local 
ULIRGs, but where the i[cii]/iFiR ratio is enhanced by 
the relatively large magnification of a localized HII re- 
gion that falls near a caustic. However, in the absence 
of other evidence supporting such a scenario, we assume 
that the measured flux ratios are equal to the intrinsic 
luminosity ratios. 

4.2. MIPS JI428 as an Extended Starburst 

The molecular gas in MIPS J1428 is exposed to similar 
UV radiation fields, and has similar excitation, as the gas 
in the central region of a typical starburst galaxy, and we 
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suggest that MIPS J1428 may be modeled as a scaled- 
up version of a starburst miclcus. As an example, we 
consider scaling up the central region of the prototypical 
starburst galaxy M82 in such a manner as to match the 
much larger luminosity of MIPS J1428, while at the same 
time conserving the strength of the characteristic FUV 
radiation field that controls the PDR emission. We as- 
sume that in a starburst environment this characteristic 
field is determined by the global star formation density, 
rather than by the purely local properties of the inter- 
actions between individual star clusters and their natal 
molecular clouds. This interpretation was adopted for 
M82 and NGC 253 based on multiple-hne PDR analy- 
sis (Wolfire et al. 1990; Carral ct al. 1994: Lord et al. 
1996), and is also motivated by the non- linear relation 
between /[cii] and /fir observed in larger samples of 
galactic nuclei (Crawford et al. 1985; Stacey et al. 1991). 

For a starburst region in which young stellar clusters 
and molecular clouds are randomly distributed, Wolfire 
et al. (1990) modeled the relationship between the aver- 
age FUV flux incident on the molecular gas (Go), and 
the size (D) and total luminosity (Xir) of the region. 
This relationship depends on the mean free path of a 
FUV photon (A), but in the limits of A < D and A > 
these models give Go oc Lir/Z?^ and Gq oc Lir/D^, re- 
spectively. Adopting Ljft. ~ 3 x 10^" Lq for the central 
D ^ 300 pc starburst region of M82 (Telesco & Harper 
1980; Joy et al. 1987), and assuming that MIPS J1428 
and M82 have the same average Go (Fig. 3), the factor of 
~1000/i~^ greater luminosity of MIPS J1428 then trans- 
lates into a size scale of D w 3/x"^/'^ — 10/i^^/^ kpc. In 
short, if the large luminosity of MIPS J1428 is produced 
by only moderate-intensity radiation, the star formation 
generating this radiation must be extended over a large 
area. 

While MIPS J1428 remains unresolved in dust con- 
tinuum and molecular gas emission, integral field spec- 
troscopy shows Ha emission extended over «0.7" («6 
kpc) (Swinbank et al. 2006). This result is consistent 
with our interpretation of MIPS J1428 as a source pow- 
ered by extended star formation, although the degree to 



which the Ha image is distorted by the foreground lens, 
or by the presence of large and spatially-varying extinc- 
tion, is unknown. Further interpretation of the connec- 
tion between the bright [CII] emission in MIPS J 1428 and 
the potentially large extent of the star formation will ben- 
efit from a size measurement in an extinction-free tracer, 
and better constraints on the lensing magnification. 

If MIPS J 1428 is indeed powered by a starburst ex- 
tending over several kiloparsecs this would be in sharp 
contrast to local ULIRGs, in which most of the emission 
traced by CO interferometry is contained within the cen- 
tral '-^l kpc (Downes & Solomon 1998). However, there 
is ample evidence that the most luminous star-forming 
galaxies at higher redshifts are more extended. The best 
studied sample of such sources are the SMGs, which like 
MIPS J1428 are submillimeter-bright galaxies with lumi- 
nosities of LiR ^ lO^'' Lq that are generated predomi- 
nantly by star formation (Alexander et al. 2005). Inter- 
ferometric imaging of the radio continuum of 2: = 1 — 3 
SMGs has shown that many of these sources have de- 
tectable emission extended over ~10 kpc (Chapman et al. 
2004) , and that the population has a median FWHM size 
of ^^5 kpc (Biggs & Ivison 2008). Similarly, many of the 
SMGs studied by Tacconi et al. (2006, 2008) have CO 
FWHM sizes of ~4 kpc. The model of extended star 
formation we suggest here for MIPS J1428 is therefore 
consistent with the large sizes of SMGs, which indicate 
that the most luminous star-forming galaxies at z = 1 — 3 
are powered by starbursts extending over much larger re- 
gions than in local systems. 
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